X 

5-S 



Classical and Quantum Complexity 
of the Sturm-Liouville Eigenvalue Problem 



2 



A. Papageorgiou 1 and H. Wozniakowski 

| 1,2 Department of Computer Science, Columbia University, New York, USA 

institute of Applied Mathematics and Mechanics, University of Warsaw, Poland 

'• February 1, 2008 

Abstract 

We study the approximation of the smallest eigenvalue of a Sturm-Liouville prob- 
■^j- \ lem in the classical and quantum settings. We consider a univariate Sturm-Liouville 

eigenvalue problem with a nonnegative function q from the class C 2 ([0, 1]) and study 
the minimal number n(e) of function evaluations or queries that are necessary to com- 
pute an ^-approximation of the smallest eigenvalue. We prove that n(e) = Q(e^ 1 ^ 2 ) in 
the (deterministic) worst case setting, and n(e) = 0(e -2 / 5 ) in the randomized setting. 
The quantum setting offers a polynomial speedup with bit queries and an exponential 
speedup with power queries. Bit queries are similar to the oracle calls used in Grover's 
algorithm appropriately extended to real valued functions. Power queries are used for 
a number of problems including phase estimation. They are obtained by considering 
the propagator of the discretized system at a number of different time moments. They 
qh allow us to use powers of the unitary matrix exp(^iM), where M is an n x n matrix 

obtained from the standard discretization of the Sturm-Liouville differential operator. 
The quantum implementation of power queries by a number of elementary quantum 
gates that is polylog in n is an open issue. 



In particular, we show how to compute an e-approximation with probability | 
using n(e) = Q(e~ 1 ^) bit queries. For power queries, we use the phase estimation 
algorithm as a basic tool and present the algorithm that solves the problem using n(e) = 
0(log £~ x ) power queries, log 2 e _1 quantum operations, and | log e _1 quantum bits. We 
also prove that the minimal number of qubits needed for this problem (regardless of 
the kind of queries used) is at least roughly ^ loge -1 . The lower bound on the number 
of quantum queries is proven in 0. 

We derive a formula that relates the Sturm-Liouville eigenvalue problem to a 
weighted integration problem. Many computational problems may be recast as this 
weighted integration problem, which allows us to solve them with a polylog number of 
power queries. Examples include Grover's search, the approximation of the Boolean 
mean, NP-complete problems, and many multivariate integration problems. In this 
paper we only provide the relationship formula. The implications are covered in |25j . 
Keywords: Eigenvalue problem, numerical approximation, quantum algorithms 
PACS numbers: 03.67.Lx, 02.60.-x 
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1 Introduction 



The study of the potential power of quantum computers has been a major theoretical chal- 
lenge. There will be an additional incentive to build a quantum computer if we can iden- 
tify computationally important problems for which quantum computation offers significant 
speedups over computation on a classical computer. 

For discrete problems, the best known quantum algorithms are due to Shor and Grover, 
see \27\ IT5] . Shor's algorithm for factorization has an exponential speedup over all known 
algorithms on a classical computer. Still, we can not yet claim that we have an exponential 
speedup for this problem, since the complexity of factorization on a classical computer is 
unknown. Grover's algorithm for data search offers a quadratic speedup. 

For continuous problems, quantum complexity is known for linear problems such as mul- 
tivariate integration, path integration and multivariate approximation, see f3J El El El 
l2l| 132] . For these problems we have an exponential speedup over the worst case setting, and 
a polynomial speedup over the randomized setting. The first quantum study of a nonlin- 
ear continuous problem was done in [20] for ordinary differential equations with polynomial 
speedups over the classical settings. 

The purpose of this paper is to present classical and quantum complexity results of 
another nonlinear continuous problem. This continuous problem is quite natural and com- 
putationally important, since it corresponds to the (simplified) univariate Sturm-Liouville 
eigenvalue problem. The Sturm-Liouville eigenvalue problem is defined in [IU] m full gener- 
ality. Here it is defined as finding the smallest eigenvalue of the differential operator 

h q u (x) := —u"(x) + q(x)u(x) for x G (0,1), 

with the boundary conditions u(0) = u(l) = 0. We assume that the function q is non- 
negative and belongs to the class C 2 ([0, 1]) of twice continuously differentiable functions 
whose norm \\q\\ := max i=0 i 2 max xe [o 1] |g^(x)| is bounded by 1. The operator h q maps 
C 2 ([0,l])intoC([0,l]). 

The Sturm-Liouville eigenvalue problem has been extensively studied in the literature. 
The properties of the eigenvalues and the eigenfunctions are well known and so are numerical 
algorithms for approximating them on a classical computer, see, e.g. [21 HE El El EH| 
Nevertheless, the complexity of approximating the smallest eigenvalue in the worst case and 
randomized settings, as well as in the quantum setting, has not yet been addressed. 

In this paper we study classical and quantum algorithms. We prove bounds on the 
worst case and randomized complexities on a classical computer, and bounds on the query 
complexity and on the qubit complexity. We prove that the complexity in the classical 
settings is a polynomial in e~ x . 

We study the quantum setting with bit queries and prove polynomial speedups over the 
classical settings. Bit queries correspond to approximate computation of function values, 
see [H] , and are used in all papers dealing with the quantum study of continuous problems. 

We also study the quantum setting with power queries. Such queries are formally defined 
in Section 5.2. Here we only mention that they are used in the phase estimation algorithm, 
which is the core of many quantum algorithms including Shor's and Grover's algorithms. 
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Power queries are controlled- W Pj queries for some n x n unitary matrix W and some expo- 
nents Pj. For the phase estimation algorithm, we have pj = 2 J_1 for j = 1, 2, . . . , m, with m 
of order log e . For the factoring problem of a large integer N, Shor's algorithm uses 
the unitary matrix W such that power queries can be implemented by at most 0(log 3 iV) 
elementary quantum gates. 

For the Sturm-Liouville eigenvalue problem, as well for all problems studied in [23], we 
use power queries with the specific unitary matrix 

PF = exp (iiM g ) with i = ^-L, (1) 

where M q is an n x n real symmetric tridiagonal matrix that is a classical approximation 
of the differential operator L g , see Section 3.2. The matrix M q depends on the values of 
q(j / {n + 1)) that appear on the diagonal of M q for j — 1, 2, . . . , n. 

Unitary matrices similar to play a key role in quantum mechanics. They give the 
solution of the Schrodinger equation, they are the propagator of a system evolving with 
Hamiltonian M q , and are important in quantum simulation, see [221 • Zalka [3^ deals with 
their implementation. The crucial point about power queries is that we can use of the 
matrix W given by as one quantum query for some j. Hence, lower bound results for bit 
queries do not apply to power queries. 

We prove that in the quantum setting with power queries, the Sturm-Liouville eigenvalue 
problem requires only roughly log e^ 1 power queries with the matrix W of (Q). As shown 
in [23], many computational problems can be reduced to the solution of the Sturm-Liouville 
eigenvalue problem, and they can be also solved in polylog number of power queries. The 
list of such problems include Grover's search, NP-complete problems, and many continuous 
problems. This proves that the quantum setting with power queries with the matrix W of 

is exponentially more powerful than the quantum setting with bit queries. 

We stress that, contrary to Shor's algorithm, we do not know if power queries with the 
nxn matrix W of (Q) can be implemented by a number of existing elementary quantum gates 
that is polylog in n. We asked a number of colleagues and most of them doubt whether this 
can be achieved. If this is indeed the case, then the positive results on the polylog number of 
such power queries will be of only theoretical interest. Still, if a future quantum computer is 
able to perform such power queries in a polylog number of, perhaps, more general elementary 
quantum gates or by some other quantum devices, the polylog number of power queries will 
lead to efficient quantum algorithms, and will allow us to solve many computational problems 
exponentially faster than on a classical computer. From this point of view, we may interpret 
the positive results on the number of power queries with the matrix W of as the indication 
that building a quantum computer with such queries would be a very desirable task which 
would give us a very powerful computational device. 

2 Survey of the Results 

In this section we explain our results in more technical terms. For a classical computer, we 
study the worst case and randomized settings in the real number model of computation with 
oracles, see j2SlE01EI]- That is, we assume that arithmetic operations (addition, subtraction, 
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multiplication, division, and evaluation of elementary functions), as well as comparisons of 
real numbers, are performed exactly with cost taken as unity. We also assume that the 
information about functions q is given by sampling q at finitely many points with the cost 
of one function evaluation taken as c. Typically c ^> 1. 

We want to approximate the smallest eigenvalue X(q) of the operator L g to within e. 
Let n(e) be the smallest number of function values of q needed to compute such an e- 
approximation in a given setting. The number n(e) is called the information complexity. The 
complexity, comp(e), is defined as the minimal total cost of computing an ^-approximation 
in a given setting. Obviously we have 

cn(e) < comp(e). 

We prove that in both classical settings, the complexity of the Sturm-Liouville eigenvalue 
problem is polynomial in e -1 , or equivalently is exponential in the number [log £^ 1 \ of 
correct bits of a computed approximation. More precisely, there exist positive numbers ol^ 
independent of e such that: 

• in the worst case setting, 

«!^ 1/2 < n(e) < a 2 e- 1/2 , 
«ice" 1/2 < comp(e) < a 2 ce~ 1/2 + a 3 e~ 1/2 log e' 1 , 

• in the randomized setting, 

a 4 e~ 2/5 < n(e) < a 5 e~ 2/5 } 
a 4 ce~ 2/5 < comp(e) < a 5 ce~ 2/5 + a e e~ 1/2 log e -1 . 

The lower bounds on n(e), and consequently on comp(e), are obtained by relating the 
eigenvalue problem to the integration problem for functions from the unit ball of C 2 ([0, 1]). 
It is well known that the minimal number of function values for this integration problem is 
bounded from below by roughly e~ x l 2 in the worst case setting and by e~ 2//5 in the randomized 
setting; see, e.g., [2311201 and the survey of these results in f3Tj . 

The upper bounds on n(e) and comp(e) in the worst case setting are obtained by the 
cost of the classical algorithm that computes an ^-approximation by the bisection algorithm 
of the Sturm sequence (SSI P- 300], see also [TTJ Ch. 5.3.4], applied to an n x n matrix which 
is the classical discretization of the operator h q with n = Q(s~ l l 2 ). The matrix depends 
on n function values of q computed at equidistant points of [0, 1]. Since we need roughly 
log £ _1 bisection steps, and the cost of each step is proportional to n, the total cost is of 
order (c + log s~ l )e~ l l 2 . Hence, modulo the logarithm of e' 1 , the worst case complexity is 
of order ce -1 ' 2 . 

The upper bounds on n(e) and comp(£:) in the randomized setting are obtained by the 
following algorithm. We first approximate the function q by a natural cubic spline q using 
n deterministic sample points of q at equidistant points of [0,1] with n = G(e -2 / 5 ). The 
relationship between the smallest eigenvalue and integration problems, see Section 3, states 
that 

X(q) = X(q) + [ (q(x)~q(x))u 2 (x)dx + 0(n~ 4 ). (2) 
Jo 
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Here Uq is the normalized eigenf unction, J Q dx = 1, corresponding to the smallest 
eigenvalue A(g). 

Since we have complete information on the spline q, we may approximate X(q) and Uq 
with arbitrarily small error. For X(q), we achieve an error of order e as in the worst case 
setting, with cost proportional to e^ 1 ! 2 log e _1 . To obtain an approximation to Ug, we apply 
one step of the inverse power algorithm with an appropriately chosen initial vector. In this 
way we obtain a vector, from which we compute Ug via piecewise interpolation. The total 
cost of computing X(q) and Uq is of order e~ 1//2 log e~ x . 

We then approximate the second term in (j2J) using the Monte Carlo algorithm for the 
function (q(x) — q{x))u^{x) computed at n randomized points with uniform distribution over 
[0, 1]. This leads to an ^-approximation in the randomized setting with cost bounded from 
above by a quantity proportional to ce~ 2 / 5 + e^ 1 ! 2 log e -1 , where the first term bounds the 
information cost and the second term bounds the combinatorial cost of the algorithm. Hence, 
we have a sharp estimate on the randomized information complexity n(e). The ratio of the 
upper to lower bounds of the randomized complexity is roughly at most 

In both classical settings, algorithms for which we obtain upper bounds on complexity 
require space of order e~ 1//2 . This follows from the fact that we need to work on n x n 
tridiagonal matrices with n of order e -1 / 2 . 

We now turn to the quantum setting. Quantum algorithms are described in Section 4. 
Here we only mention that quantum algorithms work on 2 U x 2 V unitary matrices, where v 
is the number of qubits. The qubit complexity is defined as the minimal number of qubits 
needed to solve a problem. Roughly speaking, the qubit complexity corresponds to the 
space complexity for a classical computer. For the foreseeable future, qubits will be a scarce 
resource. That is why the qubit complexity is especially important, and computationally 
important problems with relatively small qubit complexity are of special interest. 

We prove that the qubit complexity, comp qub (e), of the Sturm- Liouville eigenvalue prob- 
lem is of order log e -1 , which is relatively modest. In this paper log denotes log 2 . More 
precisely, we prove that 

\ log e~ l + < comp qub (e) < § log E" 1 + 0(1). 

These bounds hold regardless of the kind of queries used. Clearly, the qubit complexity 
yields a lower bound for the cost of any quantum algorithm solving this problem. 

We now turn to the quantum setting with bit queries. We show that the bit query 
complexity is Q(e~ 1 ^ 3 ). This result is obtained by using: 

• equation (J2J) relating the Sturm-Liouville eigenvalue problem to integration, 

• a lower bound on bit queries for integration, and 

• a modification of the classical randomized algorithm described above that uses a quan- 
tum summation algorithm instead of Monte Carlo to approximate the weighted integral 
in ©. 

We now discuss the quantum setting with power queries. In this setting, the Sturm- 
Liouville eigenvalue problem can be solved using the well-known phase estimation algorithm 
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as a basic tool, see, e.g., [221 Section 5.2]. This algorithm uses power queries and the 
quantum inverse Fourier transform as its main ingredients. The power queries have the form 
controlled- W 23 for j £ N, i.e., they use powers of the matrix W = exp (±i M q ) , with M q an 
n x n real symmetric tridiagonal matrix whose diagonal elements depend on the values of q. 
The matrix M q is a well-known discretization of the differential operator h q , and its size n 
depends on the necessary accuracy. To obtain an ^-approximation we use n of order e~ x l 2 . 

The phase estimation algorithm uses the exact eigenvector of M q , equivalently of W, as 
part of its initial state, see [221 Section 5.2]. Abrams and Lloyd pQ analyzed the case when 
the exact eigenvector is replaced by an approximate eigenvector and concluded that as long 
as the approximation is good enough, the phase estimation algorithm will still supply a good 
approximation to the corresponding eigenvalue. Jaksch and Papageorgiou ^H] proposed an 
efficient construction of an approximate eigenvector. Their idea was to solve the problem 
with low accuracy on a classical computer and obtain a "short" vector which approximates 
the eigenfunction u q at few points. Then the amplitudes of this short vector are replicated 
on a quantum computer by the Hadamard transform, which yields a "long" (vector) state 
that can be used as the approximate initial state in the phase estimation algorithm. 

We show how the construction of Jaksch and Papageorgiou can be used for the Sturm- 
Liouville eigenvalue problem. In this way, we compute an ^-approximation of the smallest 
eigenvalue with probability | by the phase estimation algorithm using log e^ 1 +0(1) power 
queries. The algorithm requires an additional number of quantum operations at most of 
order log 2 e _1 . This additional cost is for the quantum inverse Fourier transform. Finally, 
the number of qubits is § log e^ 1 + 0(1). A lower bound on the number of power queries 
of order log e~ x has been proven in [S]. 

Comparing these quantum estimates to the classical complexity bounds in the worst 
case and randomized setting, we see that the quantum setting with power queries yields an 
exponential speedup between the number of power queries and the number of function values 
needed for the Sturm-Liouville eigenvalue problem. 

Finally, we point out important consequences of our results, which we study in detail in 
[23] . Knowing that the Sturm-Liouville eigenvalue problem can be solved with polylog power 
queries, it is natural to study which computational problems can be reduced to this problem. 
In this respect, we think that the most important result of this paper is the formula that 
relates this eigenvalue problem to integration. In a particular case, this formula, see (JTUJ), 
states that 

X(q) = vr 2 + \ + 2 / (q(x) - Q sin 2 (vrx) dx + O (\\q - IQ . (3) 
Jo 

Hence, the problem of computing the smallest eigenvalue is equivalent, modulo the second 
order term, to the weighted integration problem. Since X(q) can be approximated with poly- 
log power queries, so can the weighted integral of q. It turns out that many computational 
problems can be formulated as an integration problem. Examples include important discrete 
problems such as Grover's search, the approximation of the Boolean mean, and NP-complete 
problems. The approximation of the Boolean mean is used as the primary tool to compute 
multivariate integrals and path integrals. Hence, all these problems can be solved by reduc- 
ing them to the Sturm-Liouville eigenvalue problem with a polylog number of power queries 
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in the quantum setting. It is well-known that Grover's search and the approximation of the 
Boolean mean require a number of bit queries polynomial in the problem size, which in our 
case is a polynomial in e^ 1 . This shows that power queries are exponentially more powerful 
than bit queries, see j2S] for details. 



3 Problem Definition 

We deal with functions from the class 

Q = { q : [0, 1] -> [0, 1] I qE C 2 (\0, 11) and \\q\\ := max max < 1 ). 

1=0,1,2 are [0,1] J 

For a function g 6 Q, we consider the Sturm-Liouville eigenvalue problem L, g u = Am for a 
non-zero u, or equivalently 

- q(x)u(x) + Aii(x) = 0, for x G (0, 1), (4) 
with the boundary conditions 

u(0) = u(l)=0. (5) 

Let A = X(q) be the smallest eigenvalue of (p£J), (JHJ. Multiplying (pQ) by it and integrating by 
parts, see j2 CHI I2H] , we conclude that the smallest eigenvalue satisfies 

. Jo[(u'(x)) 2 + q(x)u 2 (x)] dx 
A (o) = mm — 5 , 6 

where ifg i s the Sobolev space of absolutely continuous 1 functions for which it' G ^([0, 1]) 
and it(0) = u(l) = 0. 

Let u q be a normalized real eigenfunction corresponding to the smallest eigenvalue. It is 
known that the eigenvalues of h q are simple, and the eigenspace corresponding to A(g) is of 
dimension one. Therefore u q is uniquely defined up to the sign. In particular, u 2 q is uniquely 
defined. Then (jHJ) states that 

■M<?) = j {^{ u 'q( x )) 2 + Q( x ) u q( x ) ) dx and ||%||l 2 := (^J u 2 q (x) dx^j = 1. 

(7) 

Observe that q G Q implies that u q G C 4 ([0, 1]). Since ||g|| < 1, and ||it 9 ||z, 2 = 1 with 
u q (0) = u q (l) = 0, then \u q (x)\ are uniformly bounded for all i G [0, 4], x G [0, 1] and q G Q, 
see e.g., [TUl p. 337]. 

The smallest eigenvalue X(q) is a non- decreasing function of q, i.e., qi(x) < q2(x) for 
x G [0, 1] implies X(qi) < A(g2)- It is known that for q = c we have 



A(c) = 7r 2 + c and m c (x) = v^2 sin( 



1 A function / is absolutely continuous if and only if it can be written as f(x) = /(0) + fjf f'(t)dt for all 
as € [0,1]. 
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This implies that for q G Q, we have A(g) G [A(0), A(l)] = [n 2 , tc 2 + 1]. 

We will need estimates of the smallest eigenvalues and their eigenfunctions for perturbed 
functions q. This is a classical problem and many such estimates can be found in the litera- 
ture, not only for the simplified Sturm-Liouville problem that we consider in this paper but 
also for more general eigenvalue problems. In our case, the problem of perturbed eigenvalues 
and eigenvectors is well-conditioned, since the differential operator h q is symmetric and the 
eigenvalues of L g are well separated. Combining results from ^UJ EH I2HJ one can obtain the 
following estimates for q, q G Q: 

\X(q)-X(q)\ < ||g — g||oo := max \q(x) - q(x)\, (8) 

\\u q -Uq\\oo = O (\\q - qWoo) , (9) 

\(q) = X(q)+[ (q(x)-q(x))u 2 (x)dx + OiWq-qWl). (10) 
Jo 

We stress that the factors in the big-0 notation are independent of q and q. 

These relations follow by elementary arguments. Indeed, (JEJ) follows from (jHJ) by taking 
u = Ug, which leads to X(q) — X(q) < \\q — q\\oo- By replacing the roles of q and q we get 
X(q) — X(q) < \\q — q\\oo, which implies (JBJ). The next relation @ can be also proved by 
a matrix approximation to the operator h q , which will be done in Section 4. Finally, (fTUjl 
follows by again taking u = u q in ( J5|l . which leads to 

X(q) < X(q) + / (q(x) - q(x)) u 2 q(x) dx 
i 



o 



= X(q) + / (q(x) - q(x)) u 2 (x) dx + / (q(x) - q(x)) (ug(x) - u 2 (x)) dx. 
Jo ' Jo 

By Q, the last term is of order \\q — qW 2 ^. Taking u = u q in the expression © defining X(q), 
we obtain 

A(g) < X(q) + / (q(x) — q(x))uq(x) dx. 
Jo 

The last two inequalities imply (jlOj) . We shall see later that the formula (fTT1|) will be very 
useful in deriving lower bounds for classical algorithms. Note that if we take q = ~, then the 
formula (fTUj) becomes @. 



4 Classical Algorithms 

In this section we consider classical algorithms, i.e., algorithms on a classical (non-quantum) 
computer. These algorithms can be either deterministic or randomized. They use informa- 
tion about the functions q from Q by computing q(ti) for some discretization points U G [0, 1]. 
Here, i — 1, 2, . . . , n q , for some n q , and the points £, can be adaptively chosen, i.e., U can be 
a function 

U = ti(ti,q(ti),...,ti-i,q(ti-i)), 

of the previously computed function values and points for % > 2. The number n q can also be 
adaptively chosen, see, e.g., [SU] for details. 
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A classical deterministic algorithm produces an approximation 

<t>(q) = 0(g(fi),---,g(*n a )) 

to the smallest eigenvalue X(q) based on finitely many values of q computed at deterministic 
points. Let n = sup q£Cl n q . We assume that n < oo. The worst case error of such a 
deterministic algorithm is given by 

e wor (0,n) = sup|A(g)-0(g)|. (11) 

A classical randomized algorithm produces an approximation to X(q) based on finitely 
many values of q computed at random points, and is of the form 

<f>u,(q) = <f>Mtiju),---,Q(t nq ^)), 
where 0^, t iu> and random variables. We assume that the mappings 

LO h-> 0^, 

u i-> n g>w 

are measurable. Let n g = E(n 9ja ,) be the expected number of values of the function q with 
respect to w . As before, we assume that n = sup g6 qn g < oo. The randomized error of 
such a randomized algorithm is given by 

e ran (0,n) = sup (E[A(g) - 0.(g)] 2 ) 1/2 . (12) 
<?eQ 

For simplicity and brevity we consider the error of randomized algorithms in the L 2 sense. 
It is straightforward to extend our results for the error of randomized algorithms defined in 
the Lp-sense with p £ [1, oo]. 

We denote the minimal number of function values needed to compute an ^-approximation 
of the Sturm-Liouville eigenvalue problem in the worst case and randomized settings by 

n wor (e) = min{n : 3 such that e wor (0,n) < e } and 
n ran ( £ ) = min{n : 3 such that e ran (0,n) < e }, 

respectively. 

4.1 Lower Bounds 

We now prove lower bounds on n WOT (e) and n ran (e). 
Theorem 4.1. 

n wor (e) = n (e- 1/2 ) , n ran (e) = SI (^ 2/5 ) . 
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Proof. Define 

F={f: /GC 2 ([0,1]), maxdl/IUII/'lUlinioo) < 1 }) , (13) 
and consider the weighted integration problem 

1(f) = ! f(x)sm 2 (7rx)dx V/gF. 
Jo 

It is well-known that any algorithm using n function values for approximating of this weighted 
integration problem has worst case error at least proportional to n~ 2 in the worst case setting, 
and to n~ 2 ' 5 in the randomized setting, see fZS\ I3U] 2 . 
For c G (0, |], consider the class 

F c = F H {fEF : H/IU < c}. (14) 

For n~ 2 much less than c, the proofs for the class F can be used to deduce the same lower 
bounds on algorithms for approximation of the weighted integration problem for the class F c . 
For / G F c define q — | + /. Then q G Q. From (j3J) we have 

\(q) = ^ + l + 2 I(f) + 0(c 2 ). 

For any algorithm using n function values of q for the Sturm-Liouville eigenvalue problem, 
define the algorithm if)(f) = \(4>(q) — tt 2 — \) for the weighted integration problem. Then ip 
uses n function values of /, and 

X(q) - 0(g) = 2 (1(f) - V>(/)) + 0(c 2 ). (15) 

Let c = n~ 3 / 2 . Then n~ 2 = o(c), and therefore the error of (ft is lower bounded by VL(n~ 2 ) in 
the worst case setting, and by Q(n~ 2 ' 5 ) in the randomized setting. Hence, the error of is 
at most e when n = ^(e^ 1 ^ 2 ) in the worst case setting, and n = f2(e _2//5 ) in the randomized 
setting. Since this holds for an arbitrary algorithm 0, the proof is complete. □ 

4.2 Upper Bounds in the Worst Case Setting 

We now discuss upper bounds on n wor (e), as well as bounds on the complexity in the worst 
case setting. The worst case cost of an algorithm using n function values is defined as 

cost wor (0) = sup (c n q + m q ) , 

where m q is the number of arithmetic operations used by the algorithm for a function q 
from Q. The worst case complexity comp wor (e) is defined as the minimal cost of an algorithm 
whose worst case error is at most e, 

comp wor (e) = min{cost wor (0) : such that e wor (0, n) < e}. 

2 Formally, these results are proved for /(/) = / f(x) dx. However, the same proofs can be applied for 
the integration problem with the weight sin 2 (7ra;) and the same lower bounds hold. 
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Obviously, comp wor (» > c?i™(£). 

We now discuss the classical algorithm for the Sturm-Liouville eigenvalue problem, see 
e.g., ^UEj, and show that it is almost optimal in the worst case setting. This algorithm 
uses n = 0(e -1 / 2 ) function values of q at the equidistant points if (n + 1) for i — 1, 2, . . . , n. 
Then the operator h q is approximated by the tridiagonal n x n matrix M q of the form 



(n + l) 



2 -1 
-1 2 



2 

-1 



-1 
2 



+ 



n+l- 



Clearly, M 9 is a symmetric and positive definite matrix. Let Xj = Xj(M q ) and Zj = Zj(M q ) 
be the eigenvalues and eigenvectors of M q , i.e., M q Zj = XjZj with 

Ai < A 2 < • • • < X n , 

where the vectors Zj are orthogonal and normalized such that 



1 n 

- E z i 



fe=i 



with Zj t k being the fcth component of Zj. Note that we use the subscript L2 in the norm of a 
vector to stress similarity to the L 2 norm of functions, and to distinguish from the Euclidean 
second norm. Clearly, ||^-||l 2 = ^H^lb- 
For q = c, it is known, see, e.g., [TT], that 



and Zj(M c 



Xj(M c ) = c + 4(n + l) 2 sin 2 
[zj,i(M c ), z jt2 (M c ), z hn (M c )] T with 



JIT 



2(n+ 1) 



In 



n + l 



1/2 



sin 



jkir 
n + l 



It is known, see, e.g., [THj, that the smallest eigenvalue Xi(M q ) of the matrix M q approx- 
imates the smallest eigenvalue X(q) of the operator h q with error of order n~ 2 , i.e., 

X{q) - Xi(M q ) = O (n- 2 ) = 0(e). 

Hence, it is enough to approximate X\(M q ) with error of order e. This can be achieved by 
using roughly log e~ l bisection steps. Each step consists of computing the n terms of the 
Sturm sequence, and this can be done in cost proportional to n. The total cost is of order 
(c +log s~ l )s~ l l 2 . For details, see [TT| EH] ■ Theorem 14 . 1 1 and the cost of this algorithm lead 
to the following bounds for the minimal number of function values and for the worst case 
complexity. 
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Theorem 4.2. 

n WOT (e) = &(e- 1/2 ), ft(c£- 1/2 ) = comp wor ( £ ) = 0(ce- 1/2 + £- 1/2 log a' 1 ). 



Remark 4.1. We now show how (JHJ) can be proven, based on the properties of the 
matrix M q . First observe that for q = 0, the eigenvalues A.,- (Mo) are well separated, since 



A j+1 (M ) - A,(M ) = 4(n + l) 



(2j + 1)tt tt 



sin — ; r- sin 



2(n + l) 2(n+l) 



2 37T 7T 2 

> 4m + 1) sin — sm — « 37r . 

- V ; 2(n+l) 2(n+l) 

For q G Q, the Hermitian matrix M 9 differs from M by the diagonal matrix diagg(i/(n+l)) 
whose elements satisfy q(i/(n+ 1)) G [0, \\q\\oo] with ||g||oo < 1- Using the known estimates 
on the perturbed eigenvalues of Hermitian matrices, see [33] , we have 



min \\j{M q )- \(M )\ < Wq^ 

i=l, 2,. ..,n 

for all j = 1, 2, . . . , n. Since the intervals [Aj(M ) — 1, Aj(M ) + 1] are disjoint, we conclude 
that 

\Xj(M q ) - Xj(M )\ < IMU < 1, 

and that 

X j+ i{M g ) - Xj(Mg) > X j+ i{M ) - Xj{M ) - 2 « 3vr 2 - 2. 

Define 

J_\ /_* 
n + i;''"' 'U + l. 
where w g is the normalized real eigenfunction corresponding to the smallest eigenvalue. Then 
ll%7,nlU 2 = 1 + o(l). We normalize u q ^ n and obtain 

1 

^q,n II ~ u ^q,n- 

\\">q,n\\L2 

As mentioned in Section 3, the eigenfunction u q is defined uniquely up to its sign. Ob- 
viously, the same is true for the eigenvector zi(M q ). We choose the signs of u q and zi(M q ) 
such that 

\\u q , n - z 1 {M q )\\ L2 < \\u q , n + zi(M q )\\ L2 . 
All the components of the vector 

r] n := M q u q , n - X(q)u gin 

are of order n~ 2 , and therefore || ?7n, || z, 2 = 0(n~ 2 ). From the a posteriori error estimate, see 
|3lfl p. 173], we conclude that 

K,„-zi(M,)|| La = 0{n- 2 ) Vg G Q 
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with the factor in the big-0 notation independent of q. Note also that 

M q Uq tn - \{q)Uq,n = MqUq, n - \{q)Uq,n + T n , 

with ||r„|| i2 = 0(\\q - qW^). Hence 

||u ?f „ - Zi(M q )\\ L2 = 0(\\q- q\\oo + n~ 2 ). 

Finally, we have 

IK,n - Uq,n\\L 2 = \\Uq >n - Zl(M q ) + Z!(M q ) - Ug jn \\ L2 = O (n~ 2 + \\ q - q || «,) . 

Letting n tend to infinity, we conclude that 

\\Uq-Uq\\ L2 = 0(\\q - qWoo). 

Since both u q and Uq satisfy (jH) for (q, A(g)) and (q, A(g)), respectively, we have 
Therefore 

\W'q-ul\\ L2 = O(||g-g||oo). 
This and the fact that u — Uq vanishes at and 1 imply 

IK - %||oo = 0(\\q - g||oo), 
as claimed. □ 

4.3 Upper Bounds in the Randomized Setting 

We now turn to the randomized setting. The cost of a randomized algorithm 0, using 
n = sup g Q E(n 9)Ct! ) < oo randomized function values, is now defined as 

1 /2 

cost ran (0) = sup (E (cn ?iU + m q ^) 2 ) 

where m q ^ is the number of arithmetic operations used by the algorithm for a function q 
from Q and a random variable uj. The randomized complexity 

comp ran (e) = min { cost ran (0) : such that e ran (0,n) < e} , 

is the minimal cost of an algorithm whose randomized error is at most e. Obviously, 
comp ran (£:) > cn ra - n (e). 

We now derive upper bounds on n ra,n (e) and comp ran (e) by presenting a randomized algo- 
rithm that depends on a number of parameters. Then we find the values of these parameters 
for which the randomized error is e. We first compute m+1 function values of q at determin- 
istic points i/m, for % = 0,1, ... ,m, and construct a cubic natural spline g cu b interpolating 
q at these points, see e.g., [7j for information about cubic splines. It is well known that this 
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can be done with cost proportional to m, and \\q — g cu b||oo = 0{m~ 2 ). The function g cu b does 
not have to be non-negative. Since q > then q = q cu ^ + c > with a constant c = 0(m~ 2 ). 
We have q G Q and \\q — qW^ = 0(mT 2 ). We apply the formula (fTUj) for the function q and 
obtain 



i 



A(g) - X(q) = / (q(x) -q{x))u\(x)dx + (m~ 4 ) . (16) 
Jo 

This suggests that we can improve the accuracy of approximating X(q) — X(q) by using the 
classical Monte Carlo algorithm applied to the first term of the right hand side of (j!6|) . 
We will need to know, at least approximately, the eigenvalue X(q) and the eigenvector Uq. 
Suppose we approximate X(q) by Xq with the worst case error 

sup|A(g)-A 3 -| < S u (17) 
and the eigenfunction Uq by Zq with the worst case error 

SUp \\Uq - Zq\\ L2 < 5 2 . (18) 
<?eQ 

Assume for a moment that Xq and Zq have been computed. For a function v, define 
f v (x) = (q(x) - q{x))v 2 (x) and I(f v ) = f Q f v (x) dx. 

The randomized algorithm <p based on the Monte Carlo with k randomized samples takes 
the form 



z q{%j,u)i 



where independent and uniformly distributed numbers from [0, 1]. Here u represents 

a random element. We have 



|A(g) - <PM\ < |A(g) - Xq\ + |J(L_) - /(/ 



+ 



k 



3=1 



0(m 



Clearly, 

\Hfu,)-I(U)\ < I \q{x)-q{x)\\u\{x)-z 2 {x)\ dx = 0{rrT 2 5 2 ). 
Jo 

Since ||/z 5 ||l 2 = 0(m~ 2 ), the well known formula for the randomized error of Monte Carlo 
yields that 



* 2 ^ 1/2 




k ^,J*9\"M k l/2 



We have obtained the bound 

e ran (0, n) = O (8 X + m- 2 5 2 + m" 2 *;" 1 / 2 + m" 4 ) 
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on the randomized error of 0. Hence, to guarantee error at most e, it is enough to take 

<*i = 6(e), m = k = Q{e~ 2/5 ) and 5 2 = 9(e 1/5 ). 

We now explain how to achieve (|17|) and (fTSj) . To get Xg approximating X(q) with error 
of order e, we approximate the operator L 5 by the matrix Mg as in the worst case setting, 
now with n = Q(e^ 1 ^ 2 ). Then A(g) — Xi(Mg) = 0(n~ 2 ) = 0(e), and we compute Xg as an 
e-approximation of Xi(Mg) as for the worst case setting. This can be done with cost of order 
e -1 / 2 function values of q, and of order e~ 1 ^ 2 log e" 1 arithmetic operations. Since the cost 
of computing one function value of q is of order 1, the total cost of computing Xg is of order 
e -V 2 log £ -i. 

To get Zg approximating Ug with error of order e 1 ^ we proceed as follows. Consider the 
eigenvector z\(Mg) of the matrix Mg, with n not yet specified. By Remark 4.1, we have 

\\ug, n - Zl {Mg)\\ L2 = 0{n~ 2 ). (19) 

We approximate the smallest eigenvalue Xi(Mq) by A, with error 5. This can be achieved 
with cost of order nlog Without loss of generality we assume that A ^ Xi(Mq). Indeed, 
we can check this condition by computing det(M ? — Xl) and if this determinant is zero we 
perturb A a little. Then the matrix 

A = (Mg - XI) 1 

is non-singular and its eigenvalues are j3j = (Xj(Mg) — A) -1 . Note that \/3i\ > and 
(3j = 0(1) for j > 2. For the jth vector ej — [0, ... ,0,1,0, ... , 0] T with 1 in the jth position, 
define 

Xj = Aej. 

We can compute Xj with cost of order n by solving the tridiagonal linear system (Mg-Xl)xj = 
Cj. Then we compute 

IIeo-JU = max \\xj\\ 2 , 

j=l,2,...,n 

and 

>y I I / >"■ I I 1 ry 

1 1 JO 1 1 2 ^JO' 

Observe that the cost of computing z is of order n 2 . 

Since {n _1 / 2 Zj(Mg)}™ =1 is orthonormal, we have \\xj\\2 = Y^e=i Pl( e ji n ~ 1 ^ 2 zi(Mg)) 2 and 
|| n^ 1 / 2 ^! (Afq) HI = 1 = Xlj=i( e j' n ~ 1 ^ 2 z i(Mg)) 2 . Hence, there exists an index j such that 

(e v n~ 1 / 2 z 1 (Mg)) 2 > n-\ 

and therefore 

Ikjolh > \\ x jh > 5~ l rT l/2 . 

We have 

{Mg-X x {M q )I)z = (Mg-XI)z + (X-X 1 (Mg))z = Tr^e^ + (A - Xi(M q ))z, 

iFiolb 
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and therefore 

||(M ? -Ai(Af ? )I)z|| 2 < 5yfii + 5. 
From [SHI P- 173], we conclude that ||n -1 / 2 Zi(Mq) — = 0(8 y/n), and 

\\ Zl {M q ) - V^z\\ L2 = 0(8y/n). (20) 

We are finally ready to define Zq by piecewise linear interpolation from the successive 
components of the vector yfnz = [zi, ■ ■ ■ , z n ] T . More precisely, for j = 0,1, ... ,n let 
tj = j j (n + 1). For t G [tj, tj + i], we set 

z 9 (t) = Zj (l - (n + l)t + j) + z*+i((n + l)t - j) 

with z = z n+ \ = 0. 

We need to estimate Uq — Zq in the L 2 norm. Observe that for t G [tj, tj + \] we have 
Ug{t) = u q (tj)(l - (n + l)t + j) + u ? (tj-+i)((n + l)t - j) + 0{n~ 2 ) 
since w 5 G Q. Therefore 

|u 5 (t) - Zq(t)\ < \ug(tj) - Zq(tj)\ + \uq{t j+1 ) - Zq(t j+1 )\ + 0(n~ 2 ). 

This yields 

n ,t 3+1 
= (Uq(t) - Zg(t)) 2 dt 

j = 

Hence, 

\\uq-Zq\\ L2 = O (\\uq tn - Vnz\\ L2 + n~ 2 ) 

Since \\uq >n — \/nz\\L 2 < ||%, n — Zi(-W$)||l 2 + ||zi(Mg) ~~ V^^lUa? we use P-^l* an d to 
see that 

||%-%||l 2 = 0(5y/n + n~ 2 ). 

For 5 = n~ 5 / 2 we obtain 

11% - z q\\l 2 = 0(rT 2 ). 

Setting n = Ofe -1 / 10 ) we obtain ()18j) with 5 2 = 0(^ 1 ^ 5 ). The cost of computing Zq is of order 
n 2 = 6(e- 1 /5). 

Theorem 14.11 and the cost of this randomized algorithm lead to the following bounds on 
the minimal number of function values and the randomized complexity. 

Theorem 4.3. 

n ran (e) = 9(e- 2/5 ), n(ce~ 2/5 ) = comp Tan (e) = 0(ce- 2/5 + e~ 1/2 \og a' 1 ). 
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5 Quantum Setting 



We now turn our attention to the quantum setting. In this setting, we are using hybrid 
algorithms that are combinations of classical algorithms using function values, as explained 
in the previous sections, and quantum algorithms which we now describe. A quantum 
algorithm applies a sequence of unitary transformations to an initial state, and the final 
state is measured, see [3J |HJ EH 122] for the details of the quantum model of computation. 
We briefly summarize this model to the extent necessary for this paper. 

The initial state \ipo) is a unit vector of the Hilbert space U v = C 2 <g> ■ ■ ■ <g> C 2 , v times, 
for some appropriately chosen integer z/, where C 2 is the two dimensional space of complex 
numbers. Obviously, the dimension of 7i v is 2 V . The number v denotes the number of qubits 
used in quantum computation. 

The final state is also a unit vector of 7i u and is obtained from the initial state |^o) 
by applying a number of unitary 2 U x 2 U matrices, i.e., 



Here, Uo, Ui, . . . , Ut are unitary matrices that do not depend on the input function q. The 
unitary matrix Qy with Y = [q(ti), . . . ,q(t n )} is called a quantum query and depends on 
n, with n < 2 U , function evaluations of q computed at some non-adaptive points ti G [0, 1]. 
The quantum query Qy is the only source of information about q. The integer T denotes 
the number of quantum queries we choose to use. 

At the end of the quantum algorithm, a measurement is applied to its final state \ip). The 
measurement produces one of M outcomes, where M < 2 V . Outcome j G {0, 1, . . . , M — 1} 
occurs with probability py(j), which depends on j and the input Y. For example, if M = 2 V 
and the final state is = ^i=ci c j\j) > with X^o 1 \ c j\ 2 = 1' then a measurement in 
the computational orthonormal basis {\j)} produces the outcome j with probability py(j) = 
\cj\ 2 . Knowing the outcome j, we compute an approximation Xy(j) of the smallest eigenvalue 
on a classical computer. 

In principle, quantum algorithms may have many measurements applied between se- 
quences of unitary transformations of the form presented above. However, any algorithm 
with many measurements and a total of T quantum queries can be simulated by a quantum 
algorithm with only one measurement at the end, for details see e.g., [Tlj . 

We stress that classical algorithms in floating or fixed point arithmetic can also be written 
in the form of (J2T]) . Indeed, all classical bit operations can be simulated by quantum com- 
putations, see e.g., [I]. Classically computed function values will correspond to bit queries 
which we discuss in Section 5.2. 

In our case, we formally use the real number model of computation. Since the Sturm- 
Liouville eigenvalue problem is well conditioned and properly normalized, we obtain practi- 
cally the same results in floating or fixed point arithmetic. More precisely, it is enough to use 
0(log e~ l ) mantissa bits, and the cost of bit operations in floating or fixed point arithmetic 
is of the same order as the cost in the real number model multiplied by a power of log e -1 . 

Hence, a hybrid algorithm may be viewed as a finite sequence of algorithms of the form 
It is also known that if we use finitely many algorithms of the form (J21)) then they can 
be written as one quantum algorithm of the form (|21j) . see 




(21) 
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That is why an arbitrary hybrid algorithm in the quantum setting is of the form (|21|). 
This is important when we want to prove lower bounds because it is enough to work with 
algorithms of the form (|21|) . For upper bounds, it seems to us more natural to distinguish 
between classical and quantum computations and charge their cost differently. The cost 
of classical computations is defined as before whereas the cost of quantum computations is 
defined as the sum of the number of quantum queries multiplied by the cost of one query, 
and the number of quantum operations besides quantum queries. It will be also important 
to indicate how many qubits are used by the quantum computations. 

We now define the error in the quantum setting. In this setting, we want to approximate 
the smallest eigenvalue X(q) with a probability p > |. For simplicity, we take p = | for 
the rest of this section. As it is common for quantum algorithms, we can achieve an e- 
approximation with probability arbitrarily close to 1 by repetition of the original quantum 
algorithm, and by taking the median as the final approximation. 

The local error of the quantum algorithm with T queries that computes Ay(j) for the 
function q G Q and the outcome j G {0, 1, . . . , M — 1} is defined by 

e(Ay,T) = min j a : p Y (j) > § j- 



j: \\(q)-\ Y {j)\<a 



This can be equivalently rewritten as 



e(Ay,T) = min max | \(q) — Ay(j) | , 

where A C {0, 1, . . . , M - 1} and n(A) = E^aMj)- 

The worst probabilistic error of a quantum algorithm A with T queries for the Sturm- 
Liouville eigenvalue problem is defined by 

e i uant (A,T) = sup|e(Ay,T): Y = [qfa), . . . , q(t n )), t 4 G[0,l], for g G Q j. (22) 



5.1 Bit Queries 

Quantum queries are important in the complexity analysis of quantum algorithms. A quan- 
tum query corresponds to a function evaluation in classical computation. By analogy with 
the complexity analysis of classical algorithms, we analyze the cost of quantum algorithms in 
terms of the number of quantum queries that are necessary to compute an ^-approximation 
with probability |. Clearly, this number is a lower bound on the quantum complexity, which 
is defined as the minimal total cost of a quantum algorithm that solves the problem. 

Different quantum queries have been studied in the literature. Probably the most com- 
monly studied query is the bit query. For a Boolean function / : {0, 1, . . . , 2 m — 1} — ► {0, 1}, 
the bit query is defined by 

Qf\j)\k) = \j)\k®f(j)). 

Here v = m + 1, \j) G H m , and \k) G Hi with © denoting the addition modulo 2. For real 
functions, such as functions q, the bit query is constructed by taking the most significant 
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bits of the function q evaluated at some points tj. More precisely, as in [T3j, the bit query 
for q has the form 

Q q \j)\k) = \j)\k®P(q(T(m), 

where the number of qubits is now v = m'+m" and \j) G Tt m i, \k) G TC m " with some functions 
(3 : [0, 1] -> {0, 1, ... , 2 m " - 1} and r : {0, 1, . . . , 2 m ' - 1} -> [0, 1]. Hence, we compute q at 
tj = t(J) G [0, 1] and then take the m" most significant bits of q{t 3 ) by /3(q(tj)), for details 
and a possible use of ancilla qubits see again [T3] . 

Using bit queries, the well known quantum algorithm of Grover J3| requires Q(N 1 ^ 2 ) 
queries for searching an unordered database of N items. Similarly, the quantum summation 
algorithm of Brassard et al. [Hj computes the mean of a Boolean function defined on the 
set of N elements with accuracy e and probability | using of order minjiV, bit queries. 
Both algorithms are optimal modulo multiplicative factors in terms of the number of bit 
queries. 

The quantum summation algorithm can be also used for the approximate computation 
of the mean of a real function / : [0, 1] -»• R with |/(a;)| < M for all x G [0, 1], see UH H3- 
More precisely, if we want to approximate 

1 N ^ 

for some Xj G [0, 1] and N, then the quantum summation algorithm QS N (f) approximates 
Sjv(/) suc h that 

|Sjv(/) - QS N (f)\ < £ with probability | (23) 

using of order min(A r , Me -1 ) bit queries, min(iV, Me -1 ) log ./V quantum operations, and 
log iV qubits. 

Bit queries have been also used for a number of continuous problems such as multivariate 
and path integration, multivariate approximation, and ordinary differential equations. Tight 
bit query complexity bounds are known for a number of such problems, see [HJ EH UZ1 

In particular, Novak [21] proved that for the integration problem f(x) dx for functions 
/ from the class F given by (J13)) . the bit query complexity is 

n bit - query (£, INT F ) = Q{e~ 1/3 ). (24) 

Here and elsewhere by the bit query complexity we understand the minimal number of bit 
queries needed to compute an ^-approximation to a given problem with probability |. In 
particular, n blt_query (e) denotes the bit query complexity of the Sturm-Liouville eigenvalue 
problem. 

Based on the result ()24j) of Novak and the relationship between the Sturm-Liouville 
eigenvalue problem with integration, we now prove the following theorem. 

Theorem 5.1. 

n bit-query^ = ^( £ ~l/3)_ 
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Proof. We first prove that the bit query complexity for the weighted integration problem 
for the class F c given by (fT4"j) is of the same order as for integration for the class F, 

n bit - query (£, INTVJ = e(e~ 1/3 ). (25) 

The upper bound follows from To prove the lower bound, we use the standard proof 

technique of reducing the integration problem to the mean Boolean summation problem for 
which a lower bound on bit queries is known. 

Assume then that we use an arbitrary quantum algorithm with k bit queries that com- 
putes an ^-approximation with probability | for the integration problem over the class F c . 
Without loss of generality we assume that k~ 2 < c. 

Consider the function h(x) = ax 3 (l — x) 3 for x G [0, 1] and h(x) = for x > 1. Here, a 
is a positive number chosen such that h G F with F given by JT3J). For j = 0, 1, . . . , iV — 1, 
with N > k, define hj(x) = N~ 2 h(N(x — j/N)). Clearly, hj G F and the support of hj is 
(J/N, (j + l)/N. Observe that ||/ij||oo < N~ 2 . Hence hj G F c . We also have j hj(x)dx = 
N~ 3 Jg 1 h(x) dx. For an arbitrary Boolean function B : {0, 1, . . . , N — 1} — ► {0, 1}, define the 
function 

N-1 

f B (x) = J2 B U) h M VxG [0,1]. 

3=0 

Then f B G F c and 

Hence, modulo the factor of order A^ -2 , the computation of the Boolean mean is reduced to 
the integration problem. Note that /b(£) = B(j)hj(t) if t G [j/iV, (j + l)/iV], and sampling 
of fs is equivalent to sampling of B. From [21] we know that Q(k~ l ) is a lower bound 
for the error of the quantum approximation of the Boolean mean, with k bit queries, and 
probability |, where N > (3k for some positive (3. Letting N = \(3k~\, we conclude that the 
corresponding lower bound on the integration problem over the class F c is Q(k~ 3 ). Hence to 
achieve the error e we must have k = Q(e^ 1 ^ 3 ), as claimed in (J25|) . 

The same proof techniques allows us to consider the classes F c ( e ) with varying c(e) , even 
with c(e) tending to zero, although not too fast. We have 

n bit - query (e,INT Fc(£) ) = e(^ 1/3 ) if limc(e) e~ 2/3 = oo. (26) 

We now turn to the Sturm-Liouville eigenvalue problem. As in the proof of Theorem 14. 11 
for / G F c with c G (0, |], we define q — |+/ and consider an arbitrary quantum algorithm 
that uses k quantum bit queries and computes an e- approximation of the smallest eigenvalue 
with probability |. Then if)(f) = \((p(q) — tt 2 — \) is a quantum algorithm for approximating 
the integration problem over the class F c . We have 

\I(f)-m\ = ||(A(g)-0(?)) + O(c 2 )| < \e + 0(c 2 ). 

Take now c = c(e) = Q(e 2/3 ~ s ) with 5 G (0, |). Then 

W)-*l>(f)\ < \e + 0(e i,3 - 2& ) = |e(l + o(l)) < e for small e. 
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Hence, the quantum error of ip with probability | is e, and ip uses bit queries. Due to (j2Hj) . 
we have /c = fi(e -1 / 3 ) which completes the proof. □ 



We now derive upper bounds on the bit query complexity n blt_query (e) and on the total 
quantum complexity comp blt_quant (£:). The total quantum complexity is defined as the mini- 
mal cost of a hybrid algorithm that solves the Sturm-Liouville eigenvalue problem with error 
at most e and probability |. The hybrid algorithm may require some classical computations 
and the use of function values and the cost of them is defined just as before. It may also 
require some quantum computations and the cost of them is defined as the sum of the num- 
ber of bit queries multiplied by the cost of one such query plus the number of additional 
quantum operations. The cost of one bit query is denoted by c bit . 

We present a hybrid algorithm, which will be a combination of the classical algorithm 
from Section 4 and the quantum summation algorithm QS N for a properly chosen N. We 
proceed as in Section 4 and use the same notation. From (fTTj). and (JTBJ), we have 

Ko) = A 5 + / (<l( x ) - q( x )) z g( x ) dx + + m- 2 5 2 + m~ 4 ) (27) 
Jo 

with Si, 62 and m to be specified later. Let 

f(x) = (q(x) - q(x))zq(x) a; 6 [0,1]. 

Observe that f(x) = 0(m~ 2 ), and f(x) depends on q(x), and q{i/m) for % = 0,1,..., m, 
which are used in the construction of q. Furthermore, we can compute f(x) by computing 
one function value q(x) and one function value of the already computed functions q and Zq 
at x. We approximate J Q f(x) dx by 

3=0 V 7 

with N — (m + l)k, where the parameters m and k will be specified later. Since / is twice 
continuously differentiate and f"(x) is uniformly bounded on the subintervals (i/m,(i + 
l)/m) for i — 0, 1, . . . , m — 1, it is easy to see that 

[ f(x)dx - S N (f) = O 
Jo 

We define N such that iV -2 is of order e. 

We now apply QS N (f) algorithm to compute an 9(e)-approximation with probability | 

to Sjv(/), or, equivalently to f(x) dx. To do it, we need to use the bit query Qf for the 
function /, although so far we assumed that we can use only bit queries Q q for the functions q 
from Q. This problem is resolved in Section 2 of where it is shown that algorithms using 
the bit query Qf can be simulated by algorithms using bit queries Q g at the expense of 
multiplying the number of bit queries by a factor of 2. 

From this and (J23|) with M = 0(m~ 2 ), we conclude that its is enough to perform of order 
minfV -1 / 2 , m~ 2 e~ l ) bit queries, min(e _1//2 , m~ 2 e _1 ) log e -1 quantum operations, and using 
of order log e~ 1//2 qubits. 
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We finally approximate X(q) by the following algorithm 

0(g) = \- + Q Sjv (/). (28) 

This algorithm differs from the randomized algorithm of Section 4 since we now apply the 
QS^v quantum algorithm instead of Monte Carlo to approximate L f(x) dx. Its error is 
clearly of the form 

e bit-quant^ Jij _ Q ^ + ^-2^ + m _4 + ^ _ 

To guarantee that this error is at most e, we take 

5i = Q(e), m = Q{e~ 1/3 ), k = Q{e~ 1/6 ) and 5 2 = 6(e 1/3 ). 

Using the cost analysis of Section 4 and the results of this section, we conclude the following 
theorem. 

Theorem 5.2. The Sturm- Liouville eigenvalue problem can be solved in the quantum setting 
with bit queries by the algorithm <p defined by A2fy) . This algorithm approximates the smallest 
eigenvalue X(q) with error at most e and probability | using of order 

• e^ 1 / 3 bit queries and function values, 

• e -1 / 3 log e' 1 quantum operations, 

• e~ l l 2 log £ _1 classical operations, 

• log e^ 1 qubits. 
Furthermore, 

^bit-query = Qfc-l/^ 

and 

ft(c b it^ 1/3 ) = comp bit ™(e) = 0((c + c bit )^ 1/3 + ^ 1/2 log e' 1 ) . 

Hence, we have a sharp bound of order e -1 / 3 on the number of bit queries whereas the 
upper bound on the total cost depends, as in the worst case and randomized settings, on 
£ -1 / 2 k)g s^ 1 , which is the cost of classical computations. 

5.2 Power Queries 

In this subsection we study power queries. We formally define them as follows. For some 
problems, a quantum algorithm can be written in the form 

IV) := uJv m U m . 1 W m ^---U 1 W 1 U \^ }. (29) 

Here Ui,...,U m denote unitary^ matrices independent of the function q, just as before, 
whereas the unitary matrices Wj are of the form controlled-W^, see j2"2l p. 178]. That 
is, Wj = W Pj for annxn unitary matrix W that depends on the input of the computational 
problem, and for some non-negative integers pj, j = 1,2, ... ,m. Without loss of generality 
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we assume that n is a power of two. Let be orthonormalized eigenvectors of W, 

W\yk) = cnk\yk) with the corresponding eigenvalue a^, where \aik\ = 1 and = e lAfe with 
A fc e [0,2n) for k = 1,2, ... ,n. For the unit vectors \x e ) = ae\0) + (3 e \l) G C 2 , i = 1, 2, . . . , r, 
the quantum query Wj is defined as 

Wj \ Xl )\x 2 ) ■ ■ ■ \x r )\y k ) = \ Xl )\ ■ ■ ■ la^a) (aj\0) + ^ p ^\l)j \x j+1 ) ■ ■ ■ \x r )\y k ). (30) 

Hence, Wj is a 2 U x 2^ unitary matrix with v = r + log n. We stress that the exponent pj 
only affects jthe power of the complex number e 17Afc . 

We call Wj a power query since they are derived from powers of W. Power queries have 
been successfully used for a number of problems, see again [22], including the phase estima- 
tion problem that will be discussed in the next section. The phase estimation algorithm, 
see |H1 is at the core of many quantum algorithms. It plays a central role in the fast 
quantum algorithms for factoring and discrete logarithms of Shor [22j. We stress that for 
Shor's algorithm, power queries can be implemented by a number of elementary quantum 
gates that is polylog in n. The phase estimation algorithm approximates an eigenvalue of a 
unitary operator W using the corresponding eigenvector, or its approximation, as part of the 
initial state. The powers of W are defined by pi = 2 i_1 . Therefore, phase estimation uses 
queries with W\ = W, W 2 = W 2 , W 3 = W 22 , . . . , W m = W 2 " 1 ' 1 . It is typically assumed, 
see |Hj, that we do not explicitly know W but we are given quantum devices that perform 
controlled-iy, controlled- W 2 , controlled-W 72 , and so on. 

For the Sturm-Liouville eigenvalue problem, as well as for problems studied in [23], we 
will use the matrix 

W = exp (ijM q ) with i = \f —1 and a positive 7, (31) 

where the nxn matrix M q was introduced in Section 3.2 as a discretization of the differential 
operator h q . The matrix W is unitary since M q is symmetric. 

For the Wj with the matrix W of (|3*T|) we modify the query definition in equation (|2ip 
and assume, as in [221 Ch. 5], that for each j the Wj is one quantum query. Accordingly, 
for algorithms that can be expressed in the form (J29j) . the number of power queries is m, 
independently of the powers pj. By analogy with (J22|) . we denote their error by e quant (A, m). 

Allowing quantum algorithms of the form (|29p with power queries, we define the power 
query complexity n power_query (£:) to be the minimal number of power queries required to 
approximate the Sturm-Liouville eigenvalue problem with error e, i.e., 

n power-query^ = min { m . 3 \ SUch that e qUant (A,m) < £>. 

The cost of one power query is denoted by c power . The total complexity, comp power_query (e), 
is the defined as the minimal cost of a hybrid algorithm in the same way as for bit queries. 

We will use the phase estimation algorithm as a basic module for approximating the 
smallest eigenvalue A(g). As shown by Abrams and Lloyd pQ, the phase estimation algorithms 
can also be used if a good approximation of the eigenvector corresponding to the smallest 
eigenvalue is known. Such an approximation is obtained by the algorithm of Jaksch and 
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Papageorgiou ^K|- Combining these algorithms, we obtain the quantum algorithm that 
computes the smallest eigenvalue with error e and probability | using 6 (log e^ 1 ) power 
queries, and 0(log e~ l ) qubits. 

For the sake of completeness, we review the phase estimation problem and algorithm, 
the results of Abrams and Lloyd and the results of Jaksch and Papageorgiou in the next 
subsections. 



5.3 Phase Estimation 

Consider W defined by (pTT|) with 7 = |, i.e., 

W = exp QiM, 

The eigenvalues of W are e lXj ( Mq >' 2 with Xj(M q ) being the eigenvalues of the n x n matrix 
M q and n is assumed to be a power of two. These eigenvalues can be written as e 271 " 1 ^, where 

Vi = vA M *) = 

are called phases. We are interested in estimating the smallest phase (pi(M q ), which belongs 
to (0,1) since Xi(M q ) e [tt 2 ,^ 2 + 1]- For convenience, we renumber and normalize the 
eigenvectors of M q , and also of W, as 

\Vj) = Vn\z j+1 {M q )), 

for j — 0, 1, . . . , n — 1. We will use {\yj)} as the orthonormal basis of the space. 

Phase estimation, see [221 Section 5.2], is a quantum algorithm that approximates the 
phase !fi(M q ). Note that to compute an ^-approximation of Xi(M q ), it is enough to com- 
pute an e/(47r)-approximation of (fi(M q ). The original phase estimation algorithm has been 
derived for the initial state |0® m )|?/o), where m is related to the accuracy and will be deter- 
mined later, and 1 2/0) = \yo(M q )) is the eigenvector of the matrix M q corresponding to the 
smallest eigenvalue X\(M q ). Abrams and Lloyd jl] showed that phase estimation can still be 
used if the eigenvector \y ) is replaced by a good approximation j^o) as the initial state. 

More precisely, expanding |^ ) in the basis of the eigenvectors \yj), the initial state takes 
the form 

n-1 
3=0 

Using m Hadamard gates, we place the first register in an equal superposition, which gives 
the state 

^11 1 n-1 

1^1) = -7= Yl ■" S \ x ^\ x ^---\ x m)J2 d j\yj)- 

xi=0x2=0 i m =0 j=0 

We now apply the controlled quantum gates, see (|30|). to create the state 



|-0 2 ) = W 2 m-lW 2 m-2 ■ ■ ■ W 2 



^ 11— 1 

==^2 d sh)\yj) 



j=0 
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with 

\ Vj ) = ( |0) + e 27ri ^|l> ) ® ( |0> + e 2 ^|l> ) ® • • • ® ( |0) + e 27ri2m "^|l) 



1 1 1 



E E • • • E e^^ 120 ^ 22 ^-^ 2 " 1 " 1 )^^!)^^ • • • |x m ) 

i'i=0 X2=0 £m=0 
2 m -l 



£=0 

see also J221 P- 222]. Hence, 

n-l /2 m -l 



Ora 

The inverse Fourier transform performed on the first register creates the state 

n-l /2 m -l \ 

where 

sin(7r(2 m y J -£))e W J . f , n-mtf 

2"»sin(7r(^-fi-"»)) 11 W r= z 

1, if </?j = 2" m £. 

A measurement of the first register produces the outcome j with probability 

n-l 

Pj = E l^l 2 b(^'i)l 2 ' 

1=0 

and the second register collapses to the state 

n-l 



The quantity 

A(0 o ,0i) = min{|x + 0i - O |} for o ,0ieR 

is defined in |H] and is the fractional part of the distance between two phases <fto and 0i. It is 
used to derive the relationship between the approximation error and the success probability. 
A measurement of the first register produces an outcome from the set 

G k = {j:A(j/2\MM q ))<k/2 m }, 

where k > 1, with probability 

= EEi^'^l 2 > \d\ 2 J2\9(Vi(M q ),j)\ 2 > \d\ 2 -J^y 

25 



where d — (yolV'o)- For k = 1 the probability that 

A(j/2 m , (pi(M„)) < T m is bounded from below by 4r|d| 2 . (32) 

The proof of the probability bounds can be found in [Ul |22]- Using this fact, the authors 
of P conclude that as long as |<i| 2 is large enough or, equivalently, |^ ) is close enough to 
\yo) then phase estimation can be used to approximate the phase ipi(M q ) with probability 
close to 8/vr 2 = 0.81.... 

We stress that the phase estimation algorithm uses m power queries. In addition to the 
cost of the queries there is a quantum operations cost proportional to at most m 2 , which is 
an upper bound on the cost of the quantum inverse Fourier transform, see |22[ Section 5.2]. 



5.4 Eigenvalue and Eigenvector Approximation 

The results of Jaksch and Papageorgiou [IE] can be applied to efficiently construct a good 
approximate eigenvector when W = e^ Mq as in the previous subsection. 

The matrix M q = M q has been derived from the discretization of the operator L g with 
mesh size h n = (n + Its eigenvectors are also eigenvectors of W — W^ n \ and we 



denote them here by {y^}, where j = 0,l,...,n — 1. We want to approximate Xi(M q n> ) 



A7C(pi(Mq n ^) but we do not know the corresponding eigenvector 



\v w ) ■= \y { * ] )- 

The expansion of \y^) in the computational basis is denoted by 



n-l 



|y (B) > = £yj n) b>, 

3=0 



Recall that u q is the normalized 



1/2 



(33) 

1, eigenfunction of the dif- 



U q\\L 2 = \yJo U 2 (x)dx^ 

ferential operator h q that corresponds to X(q), and u q as well as u' q and u' q are uniformly 
bounded, i.e., l|w,j||oo, H^Hoo and lltil'lloo are 0(1). 



Let {U^} = Y^i=o u q{(j + be the vector obtained by sampling u q at the dis- 



cretization points. Then it is known, see f21E2] as well as Remark 4.1, that 

,c»h \ uin) ) 



\v w ) 



Ill^ (n) )ll2 

\m - \x{M<?>)\ 



o{K) 
o{K) 



and 



(34) 



Consider a coarse discretization of L 9 with mesh size h no 
of two. Assume that 



(n + 1) 1 with n being a power 



no-l 

E 

3=0 



Z 3 



li>, 



|5 ( " 0) )|| 2 
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approximates the eigenvector \y^ n °') that corresponds to the smallest eigenvalue of the matrix 
il4 no) such that, 

|| |£(no)} _ \ y (n )j || 2 = O(n 2 ). (35) 

We place the vector \z <yn ° S} ) in a log Uq qubit register. As explained in Section 4.3, we can 
compute \z <yn ° S) ) on a classical computer with cost of order n 2 ,. 

For n = 2 s uq, we construct an approximation \z^) of \y^) by first appending s qubits, 
all in the state |0), to and then performing the Hadamard transformation on each one 

of these s qubits, i.e., 

ISM^IiMjp^-^g^, (36 ) 

where z^)'s denote the coordinates of {z^} in the computational basis, and g(j) = [j/2 s J. 

The effect of g is to replicate 2 s times the coordinates of \z <Jl °' > ). As in Jaksch and Papageor- 
giou [18], we use the vector \z^ n ') as part of the input to the phase estimation algorithm. 

Let d< n ) = (y^\z^). We show that \d^\ 2 can be made arbitrarily close to one by 
choosing a sufficiently large n . Hence, we can make the success probability of the phase 
estimation algorithm at least equal to |. 

Consider two different expansions of | 

n-l 

= ^ fi W U> (37) 
i=o 

n-l 

| 5 (n) } = Y^df^yf). (38) 
i=o 

The first expansion is in the computational basis {\j)} and, by (jHEJ), 

^ =2 -s/2 z M forj = ,l,...,n-l, 

while the second expansion is with respect to the eigenvectors of Mq . Note that = d^ 
and clearly Y^jZo Mj™^ 2 = 1- Equation (fH£J) implies 

n-l 

- | y W> = (rfW - l)| y W> + y £df ) \yf l) ). (39) 

i=i 

Taking norms on both sides we obtain 

n—l n— 1 

|||y(»>)-|zW)||* = |rfW-l| 2 + ^|4")| 2 > ^|4 n) | 2 = l-\d^\ 2 . (40) 
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We now bound the left hand side of (|40|) from above. Using the expression ()33|) for \y^) 
and the definition of \z^ n >), see (J31)]) . (|3*7jl . we have 



n-1 

_ 2 - s /2 z («o)|2 

j=0 
n-1 

E 

J=0 



u q {{j + l)h n ) u q {{g{j) + l)h 



no J 



\\UM)\\< 



^]||C/M)| 



A 



(n) 



ah) 



where, by (J3*4*j) and (J33J), we have 



n-1 

i=o 



n-1 



(«)|2 



O(^) and 



£|A^| 2 = 2*0(/, 4 



i=o 



Applying the triangle inequality, we get 



'n-1 



| y W>-|zM>|| 2 < [J2 
0=0 



u g(U + 1)M u q ((g(j) + l)h 



no 1 



^|||C/M)| 



no/' 



1/2 



»()' 



(41) 



The definition of If/ 1 -™)) and the fact that the derivative of u q is Lipschitz 3 with the uniform 
Lipschitz constant imply that |||C/^)||2 = + 0(h n )). Hence, the square of the term in 

the parentheses above is equal to 

1 n— 1 

- Y\ MU + + 0{h n )) - u q ((g(j) + l)h no )(l + O(h no ))\ 2 . (42) 



3=0 

Since u q is continuous with a bounded first derivative, we have that 

u q (x 2 ,j) = u q (xu) + 0(\x 2 j - x 1J \), (43) 

where x 2>j = (j + l)h n and x hj = (g(j) + l)h no , j = 0, 1, . . . , n - 1. Let [j/2 s \ = j/2 s - a 
with a G (0, 1). Then 

j + 1 3 /2 s + 1 - a 



\%2,j — I 



J 



2 s n + 1 n + 1 

2 s - 1 



Using 



(2*n + l)2«(no + 1) 
and the triangle inequality, we obtain from ()41|) that 



y in) ) - \z [n) ) 



n + 1 

for some positive number c independent of n and no- Combining this with (f4T)j) we finally 
conclude that 



n)|2 



> 1 



K + l) 2 

and S-^ can be made arbitrarily close to one by taking a sufficiently large rig. 



(44) 



3 A function / : [0, 1] 
all x,y £ [0, 1]. 



is Lipschitz if there is a number L > such that \f(x) — f(y)\ < L\x — y\ for 
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5.5 Quantum Algorithm for the Smallest Eigenvalue 

We combine the results of the previous two subsections to derive a quantum algorithm for 
computing an e- approximation of the smallest eigenvalue with probability |. 

We choose the parameters for the phase estimation algorithm. Without loss of generality 
we assume that e" 1 is an even power of 2, that is e~ l = 2 m with an even m. We set 
n = e~ 1 / 2 = 1 m l 2 and we will be working with the matrix Mq n \ The index n = 2 k ° is 
chosen as the smallest power of two for which 

Hi t 1 - (£+1)0 -4' (45) 

where the number c is from (jUJ). Clearly, no = 0(1). Without loss of generality we assume 
that \m > ko = log no, i.e., we assume that e is sufficiently small. We finally set s = ^m — k^. 
We then compute \z <yn °^ > ) on a classical computer as in Section 4.3 with cost 0(1) function 
values and operations. 

i ("■) 

We run the phase estimation algorithm for the matrix W = e? Mq with the initial state, 
see (JSnj), 

|0)®"*|g(*)) = | )® m |^ n °)) (M-^M^j . 
Let j be the outcome of the phase estimation algorithm. We finally compute 

Xj = 4tt j 2- rn 

as an approximation of the smallest eigenvalue X(q). We have 

Xj - X(q) = Xj - X l (M^) + X l (M^) - X(q) 

= ^ (^-^( M « (n) )) 

From (}3*2*|) we know that 

- ^(Mf')J < e with probability 1 \d {n \ 2 . 
By (|44|) and the definition of n Q we have 

-ilrfwi 2 >lfi- c2 1 > 2 

7T 2 ' ' " 7T 2 V K + l) 2 / 4' 

Hence, 

3 

\Xj — X(q)\ = 0(e) with probability at least -. 
The computation of Xj requires 

m + k + s = | m = | log s~ l 
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qubits, m = log e 1 power queries, plus a number of quantum operations proportional to 

m 2 = \og 2 e-\ This yields n P°wer-query( £ ) = ( log £ -ly A lower bound On n P°wer-query( £ ) Q f 

the same order is proved in [Sj. Hence, 

n power-query^) = Q ^ Qg £ -ly 

We summarize the results of this section in the following theorem. 

Theorem 5.3. The Sturm- Liouville eigenvalue problem can be solved in the quantum setting 
with power queries by the phase estimation algorithm applied to the discretized matrix of the 
differential operator L g with the initial state given as an approximate eigenvector computed by 
the Jaksch and Papageorgiou algorithm. This quantum algorithm approximates the smallest 
eigenvalue X(q) with error e and probability | using 

• | log e^ 1 + 0(1) power queries, 

• 0(1) function values and classical operations, 

• 0(log 2 e _1 ) quantum operations besides the power queries, and 

• |loge _1 +0(1) qubits. 
Furthermore, 

n power-query = QQ Qg £ -l^ 

and 

fi(c powcr log e~ l ) = comp power -™(e) = O (c power log e~ l + c + logV 1 ) . 
5.6 Qubit Complexity 

In this section we address the qubit complexity, comp qub (e), which is defined as the minimal 
number of qubits required to approximate the smallest eigenvalue with error e and prob- 
ability | by quantum algorithms of the form (J21|) . Clearly, comp qub (e) is upper bounded 
by flog e- l + 0(l) since that many qubits are used by the phase estimation algorithm of 
Section 5.5. Observe that the cost of the classical algorithm computing \z <yn ° S} ) as well as its 
quantum simulation [221 P- 189-193] is constant since no is bounded by a constant due to 
© . 

We turn to a lower bound on comp qub (e). Based on the results obtained in this paper, it 
is easy to see that the number of qubits necessary to solve our problem must be proportional 
at least to roughly |loge _1 . Indeed, assume that there is a quantum algorithm of the 
form (PTj) that computes X(q) with error e and probability |, and uses k(e) qubits. This 
algorithm can use arbitrary quantum queries, assuming that each quantum query is based on 
at most 2 fc ( £ ) function evaluations of q. Note that this holds for bit queries, as well as for the 
power queries studied in this paper. Then such an algorithm can be simulated by a classical 
algorithm that uses at most 2 k ^ function evaluations of q. From Theorem 3.2 we know that 
2k(e) _ ^(^-V 2 ) anc l therefore k(e) > |log e^ 1 + fi(l). Hence, the qubit complexity is lower 
bounded by |log e^ 1 + This proves the following theorem. 
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Theorem 5.4. The qubit complexity of the Sturm- Liouville eigenvalue problem in the quan- 
tum setting with bit or power queries is bounded by 

| log E" 1 + 0(1) < comp qub (e) < § log e^ 1 + 0(1). 
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